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METHOD FOR DETERMBNT Na STr:Nmr ANT BONE DE NSITY LOSSES 

The invention relates to a method for determining significant bone 

density losses. 

In the interaction of bone destruction processes and bone formation 
processes, there may, on the average, be loss of bone substance and structure. In 
order to determine this loss and identify significant deviations from age-specific and 
sex-specific reference values, bone density measurements (osteodensitometry) are 
generally carried out. Using the methods of photon adsorption or computer 
tomography, the bone mineral density is measured at the radius, at the thigh bone 
(femur) or at the lumbar vertebral column (Vertebra lumbalis). Under the assumption 
of an exponential development of the bone density decrease, the loss can be estimated 
from bone density values measured at least at three different times. Intra-individual 
measurement fluctuations of up to 5 percent may be obtained. In order to keep these 
factors, which interfere with the determination of the results, as small as possible, 
sufficiently large intervals between the times of measuring the bone density should be 
selected, so that changes in the bone density can be detected adequately. At the 
earliest, therefore, information is available after one year. 

Osteodensitometric methods of identifying bone density losses are 
expensive, cost intensive and expose the patient to radiation. Moreover, 

osteodensitometry is not transferable between different instruments and measurement 
methods are not standardized. At the present time, there are no recognized alternative 
methods with equivalent informative capability. 

In the literature (P. Gamero and P.D. Dalmas, "Biochemical Markers of 
Bone Turnover" 1988, Endocrinology and Metabolism Clinics of North America, vol. 
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27, No. 2, pp. 303 -322), laboratory parameters, which are particularly associated with 
bone density losses (bone markers), are investigated extensively with respect to their 
suitabiUty for describing bone density loss quantitatively. Such parameters from the 
cell portion (osteoblasts), from the organic matrix of the bone (collagen, 
noncollageous proteins) or from the inorganic constituents of the bone admittedly can 
give an indication of bone loss. Until now, however, the information, provided by 
them, is not equivalent to the information provided by osteodensitometry. With that, 
the advantage of simplicity of this method of determining parameters in the serum or 
the urine cannot be utilized. 

It is an object of the invention to develop a method for determining 
significant bone density losses, which is less cost intensive, does not expose the 
patient to radiation and shortens tiie period for interpreting the interaction of bone 
destruction processes and bone formation processes. 

Pursuant to the invention, this objective is accomplished owing to the 
fact that measurement values of real or mathematically simulated bone density loss 
processes, which are present in electronic storage media and, as a function of time, 
reflect laboratory parameters of practically or theoretically known clinical signs and 
symptoms, are used as reference values for the process, that bone marker values of 
serum or urine samples, associated with bone density losses and measured by 
common laboratory techniques, are determined by sample preparation steps such as 

• treating with antibodies 

• incubation steps 

• separation procedures 

• using analytical techniques 

and recorded over an input mask in an electronic data memory, are used for 
determining significant bone density losses, 

a) at the time of the analysis, all N available patient-related data being copied 
from the data memory over an interrogation ftmction and made available for 



the further processing (measured values M (t„; k) of the K in the laboratory of 
the bone marker, determined after step x of the process at times ti . . . tn) 

the measured values of the bone markers with respect to the first line in the 
Table being normalized according to the equation 

and the measurements as a function of time being converted into months, 

the normalized measured value being converted into a scalar quantity D(tn) for 
the graduated description of the course of the bone density, the equation 

being used as functions for the graduated description of the course of the 
relationship 

from the evaluations of the progress determined, evaluations of the progress 
for those time sections of 

^•(,)== <'--"-^<°-";<'-'--> °('-> ,,.[,..„,.] 

being calculated by interpolation, for which reference values are available, 



from the interpolated evaluations of the progress, similarity dimensions being 
calculated, the function 



M , 
11=1 ^A/ 



being used to calculate a similarity dimension between the data, which is to be 
investigated, and all the reference values, available in the database and, at the 
same time, similarity dimensions to the reference values and to the time in 
months being found; 

■L 

from the similarity dimensions for all reference values, those reference values 
being determined, which have a high similarity in the mathematical sense, 
such as tihie 



greatest similarity: A* = min{ Aj} 



positive alternative (+) a*= min {Aj} 



negative alternative (-) A"= min [Aj] 



with subsequent output of the type description as text component for 
describing the situation; 



the prediction being derived from these three reference values, if Bi = A*, B2 
= A^B3 = A-, 



the following ejcpression 




being used for the predicted value at time t 

the degrees of freedom for the specification of the model, given as functional 
parameters in the functional relation of D(t„) and Aj(t) being occupied by 
standard specifications and fitted by statistical analysis of the reference values 
to practical experience for optimizing the quantitative prediction of the bone 
density loss; 

the time being calculated, at which, according to this prediction strategy, the 
percentage deviation is greater than a specided threshold value, this time being 
the starting point for planning the scheduling of the next investigation. 

Advantageously, the degrees of freedom, given as function parameters in the 
functional relationship of D(t^ and Aj(t), are filled in by the mathematical 
method of least squares so that specified sequences are taken into 
consideration in the best way for reference values. 

The reference values used can be values from an assumed analytical 
mathematical course (exponential function), from empirically based values 
from imaginary, assumed processes and from concrete, measured values from 
patients with known situations. 



In the following, the invention is explained in greater detail by means of an 
example. 

Osteocalcin, parathyroid hormone and alkaline phosphatase are used as bone 
markers. Data for the bone markers are obtained from serum or urine samples 
using well-known laboratory techniques (HPLC, RIA, ELISA). For this 
purpose, sample preparation steps such as 

- treatment with antibodies 

- incubation steps 

- separation methods 

- use in analytical techniques 

are necessary in order to obtain, after the measurement, a quantitative value in 
the parameter-specific unit as a monitor or printer output or as an 
electronically available numerical value. These laboratory values obtained are 
stored in an electronic database over an input mask. It is a prerequisite of tiie 
method that reference values are known. Reference values can be values 
calculated theoretically from an analytically mathematically assumed course 
(exponential function) or empirically based values from imaginary, assumed 
processes or concrete values from patients with known situations. These 
reference values are available for specified times and can be used in the 
analysis only within the scope of the time horizon determined therewith. In 
the sample, the exponential function R (t)= a(l- e"^^) is used to describe bone 
density losses. The parameters a and b then describe the speed and the degree 
of loss (t in month). Some reference examples are listed in the following 
Table: 
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Months After 


Type 0 


Type I 


Type II 


Type 


Type 


Type V 


Type 


First 


a = 0 


a=l 


a = 0.5 


m 


IV 


a=l 


VI 


Measurement 




b = 


b = 


a =0.5 


a = 2 


b = 


a = 2 






0.01 


0.05 


b = 0.1 


b = 
0.01 


0.05 


b = 
0.05 


0 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


6 


0.00 


0.06 


0.13 


0.23 


0.12 


0.26 


0.52 


12 


0.00 


0.11 


0.23 


0.35 


0.23 


0.45 


0.90 


18 


0.00 


0.16 


0.30 


0.42 


0.33 


0.59 


1.19 


24 


0.00 


0.21 


0.35 


0.45 


0.43 


0.70 


1.40 


30 


0.00 


0.26 


0.39 


0.48 


0.52 


0.78 


1.55 


36 


0.00 


0.30 


0.42 


0.49 


0.60 


0.83 


1.67 


42 


0.00 


0.34 


0.44 


0.49 


0.69 


0.88 


1.76 


48 


0.00 


0.38 


0.45 


0.50 


0.76 


0.91 


1.82 


54 


0.00 


0.42 


0.47 


0.50 


0.83 


0.93 


1.87 


60 


0.00 


0.45 


0.48 


0.50 


0.90 


0.95 


1.90 



This is followed by the steps below 
a) the values measured and those obtained on three prior occasions are given 
in the Table below: 
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No. 


Date 


Osteocalcine in 


PTH 


AP 






|xg/L 


inng/L 


inU/L 


1 

1 


5-30-96 


9.8 


24.8 


90 


2 


1-29-97 


10.9 


34.6 


86 


3 


2-16-98 


12.6 


32.0 


104 


4 


3-2-99 


12.4 


34.0 


107 



measured values M(t„;k) for n = 1, . . .,4 and k= 1, . . 3 



b) the values obtained by measurement are normalized with respect to the first 
line in the Table according to the following equation 



and the time interval between the measurements is converted into months; 



No. 


Month 


M*(t;l) 


M*(t;2)^ 


M*(t;3) 


1 


0.0 


0.00 


0.00 


0.00 


2 


8.0 


0.11 


0.40 


-0.04 


3 


20.6 


0.29 


0.29 


0.16 


4 


33.1 


0.27 


0.37 


0.19 



Normalized values measured M(tn;k) for n = 1, . . .,4 and k?= 1, . . 3 
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the normalized measured value is converted into a scalar quantity for the 
graduated description of the bone density loss, the equation 



being used (K = 3; n = 1,. . .,3) as a function of the graduated description of the 
progress. Under the standard setup, it is, of course assumed that w = 1 for all 
weighting factors; 
Evaluations of Progress at Time t 

From the evaluations of the progress obtained, evaluations of the progress for 
those time intervals of 

^ ; ; '^elVij^Hj 

are calculated by interpolation, for which reference values are available; 

Interpolated Evaluations of Progress at Imaginary Measurement Times 
at 6-Month Intervals. 

from the interpolated evaluations of the progress, similarity dimensions are 
calculated, the function 



being used for calculating a similarity dimension between the data, which is to 
be investigated, and all reference values available in the database, the 
following sraiilarity dimensions being found. Under a standard setup, V = 1 
for all weighting factors. 




M 



Month 




TypeO 


Type I 


Type II 


Typem 


Type IV 


Type V 


Type VI 


0 


















6 


0.32 


0.10 


0.07 


0.04 


0.01 


0.04 


0.00 


0.04 


12 


0.43 


0.24 


0.14 


0.06 


0.01 


0.06 


0.00 


0.24 


18 


0.44 


0.36 


0.17 


0.06 


0.01 


0.05 


0.02 


0.72 


24 


0.47 


0.48 


0.19 


0.06 


0.01 


0.04 


0.07 


1.41 


30 


0.49 


0.63 


0.21 


0.06 


0.00 


0.03 


0.14 


2.25 



Similarity Dimensions for all Reference Values and all Times up to 30 Months 
The instructions for the calculations explicitly do not represent a function for 
estimating the parameters a and b of the exponential function, since general 
reference values do not have to satisfy these instructions; 

f) from the similarity dimensions for all reference values, those reference values 
are determined, which have a high similarity in the mathematical sense, such 
as the 



greatest similarity: 



= 0,00 



positive altemative (+) mill {Aj} =0,03 



negative altemative (-) A*= mm {Aj} =0,06 



In accordance with these selection instructions, Type III is recognized as the 
reference value and Type IV and Type II are selected as positive and negative, 
altemative reference value respectively; 
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g) The prediction is derived from these three reference value changes. If, in order 
to simplify the notation, 

Bi = A*, B2 = A^, B3 = A\ the following is used as the predicted value at time 
t 




1^1 



h) The degrees of freedom for the specification of the model, given as functional 
parameters in the functional relation of D(tn) and Aj(t), are filled in by 
reference values, in order to achieve a quantitative prediction of the bond 
density loss; 



Month 


Prediction for 


Prediction 


Prediction 


Predicted 


Deviation 




Typem 


for 


for 


Value R(t) 


jfrom 






Alternatives 


Alternatives 




Prediction 






(+) 


(-) 




in% 


36 


0.49 


0.60 


0.42 


0.52 


6.1 


42 


0.49 


0.69 


0.44 


0.55 


12.2 


48 


0.50 


0.76 


0.45 


0.58 


16.0 


54 


0.50 


0.83 


0.47 


0.61 


22.0 


60 


0.50 


0.90 


0.48 


0.63 


26.0 



Predicted Value 

The certainty of this prediction is characterized by the percentage deviation of 
the average predicted value R(t), using only the suitable reference value, 

ascertained with the similarity dimension factors Ai, in the example, 
prediction for 

Type ni, which has the greatest similarity A* 
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i) The time is calculated, at which, according to this prediction strategy, the 
percentage deviation is greater than a specified threshold value, this time being 
the starting point for planning the scheduling for the subsequent investigation. 
For this purpose, this value is put out ova: a monitor or printer and transmitted 
by remote data transmission to the treating physician. In the case of the 
present example, the threshold value is assumed with a 10% deviation. Forty 
months after the first measurement, a deviation greater than 10% is then found 
for the first time. The time for the next measurement should therefore not be 
later than 7 months after the values for the bone markers are determined. 



Month 


Prediction for 


Prediction 


Prediction 


Predicted 


Deviation 




Typeffl 


for 


for 


Value R(t) 


from 






Altematives 


Altematives 




Prediction 






(+) 


(-) 




in% 


40 


0.49 


0.66 


0.43 


0.54 


10.0 



Deviation of the predicted value greater than 10%. 
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